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ABSTRACT 

In this paper we present a detailed numerical investigation of the hypothesis that a rotation of 
astrophysical jets can be caused by magnetohydrodynamic shocks in a helical magnetic field. Shock 
compression of the helical magnetic field results in a toroidal Lorentz force component which will 
accelerate the jet material in toroidal direction. This process transforms magnetic angular momentum 
(magnetic stress) carried along the jet into kinetic angular momentum (rotation). The mechanism 
proposed here only works in a helical magnetic field configuration. We demonstrate the feasibility 
of this mechanism by axisymmetric MHD simulations in 1.5D and 2.5D using the PLUTO code. In 
our setup the jet is injected into the ambient gas with zero kinetic angular momentum (no rotation). 
Different dynamical parameters for jet propagation are applied such as the jet internal Alfvcn Mach 
number and fast magnetosonic Mach number, the density contrast of jet to ambient medium, or the 
external sonic Mach number of the jet. The mechanism we suggest should work for a variety of jet 
applications, e.g. protostellar or extragalactic jets, and internal jet shocks (jet knots) or external 
shocks between the jet and ambient gas (entrainment). For typical parameter values for protostellar 
jets, the numerically derived rotation feature looks consistent with the observations, i.e. rotational 
velocities of 0.1-1% of the jet bulk velocity. 

Subject headings: shock waves ~ MHD - ISM: jets and outflows - stars: mass loss - stars: pre-main 
sequence galaxies: jets 



1. INTRODUCTION 

Astrophysical jets are highly coUimated beams of high 
velocity material, observed in a variety of astronomi- 
cal sources - among them young stellar objects (YSO), 
micro-quasars, or active galactic nuclei (AGN). Some- 
what less coUimated beams of comparatively lower speed 
are usually called outflows. The current understanding of 
jet formation is that these outflows are launched by mag- 
netohydrodynamic (MHD) processes in the close vicinity 
of the central object - an acc retion disk surrounding a 
protostar or a com pact object ( Blandford &: Payne 1982} 
■Pudritz & Norman T983{ |(Jamenzind|199U[|Pudritz et al. 
2007j |Cabrit„2007^ . 

Over the past tew years observational indication has 
been accumulated for axially symmetric differences in 
the radial velocity profile across protostellar jets and out- 
flows. These radial velocity differences have repeatedly 
been interpreted as outflow rotatioi^ For one of the best 
studied examples, the micro-jet from DGTau, the sug- 
gested rotational velocity is about 6-15 km/s at about 
0.2 arcsec or 40 AU fr om the jet axis up to distances of 
100 AU from the star ( [Bacciotti et al.||2002[ ). 

As MHD outflows are launched trom rotating sources 
by nature, the rotation of the outflow itself seems to be 
a characteristic feature of MHD jet format ion models 



from classic steady s tate models ([B landford fc Payne 
1982), non-relativistic ( |Ouyed fc PuTritz|jl997| ICasse fc: 



Keppens 2002) or recent relativistic simulations! Porth 
fc I'encTEIpOlUl . The transverse stability of relativistic 



[fendt@mp ia.de| 

^ in tile following, for simplicity, we will omit such cautious 
wording regarding the existence of rotation in jets/outflows - be- 
ing aware, however, that outflow rotation has not yet been fully 
confirmed 



jets has been treated by two-component simulations by 



Meliani fc Keppens| f2007 ). When jet s originate in the in- 
ner part of a rapidly rotating disk, the outflow material is 
launched with a substantial amount of angular momen- 
tum. Magnetohydrodynamic models further imply an- 
gular momentum conservation along the streamlines of 
the outflow, and, thus, jet material which is persistently 
rotating. Quantitative estimates of launching-induced 
jet rotation based on steady-state MHD models seem to 
be in agreement with the observations ( jAnderson et al. 



20031 



Tnere are several arguments why this picture may not 
hold in general. The main argument is that the radial ve- 
locity profile is usually inferred from emission lines from 
shocked gas. It is therefore questionable that classical 
steady-state MHD models can be used to fit the veloc- 
ity and angular momentum balance along the jet across 
these shocks. Furthermore, traditional non self-similar 
steady-state models of MHD jets indicate on both large 
initial jet opening angles and large jet radii compared 
to the Alfven radius ta. In that the rotational 

velocity decreases with ~ 1/r for r > ta, jet radii of 
10^ — lO'^rA would result in a negligible rotational veloc- 
ity along the outer layers of a jet. 

We like to note here that the outfiow velocity profile 
itself might be explained without implementing jet rota- 
tion at all, and, thus, without assuming a M HD launch- 



ing model for the outflow. (Soker 



2005) suggested a 



model in which a thick and twisted accretion disk in- 
teracts with the outflow by entrainment. Inclination be- 
tween the disk rotational axis and the jet axis would lead 
to a different amount of deceleration on both sides of the 
jet, thus visible in a velocity gradient across the jet. 

In this paper, we propose a different explanation for the 
observationally indicated outflow rotation. The model 
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combines two basic features of MHD outflows - (i) the 
fact that the observed emission is from shocked gas, and 
(ii) the jet shock compressed hehcal field leads to a mag- 
netic torque which accelerates the outflow material in 
toroidal direction. Essentially, the jet kinetic angular 
momentum is retrieved from the magnetic angular mo- 
mentum content of the jet across the shock. First hints of 
such a mec hanism have been o bser ved in numerical sim - 
ulations by iKossl et al.] ( |1990[ ) or |Uchida et al.| ( |1992[) ; 



To do et al. 



1992 ), but were not discussed in great detail 
or even applied to observationt]^ 

In the following, we first summarize the observational 
findings concerning outflow rotation in jets from young 
stars and then briefly mention the case of extragalactic 
jets. We continue by discussing the governing MHD pro- 
cesses involved in our model, and present results of our 
MHD simulations. These simulations are preliminary in 
the sense that our emphasis is on the dynamical evolution 
of the outflow. This did not yet include further physi- 
cal processes such as cooling (apart from one example 
simulation), magnetic diffusivity, or radiative effects. 

2. OBSERVATIONAL INDICATION OF OUTFLOW 
ROTATION IN YOUNG STARS 

Here we summarize the observational background dis- 
cussing a few typical, but different example sources where 
outflow rotation is observationally indicated. 

2.1. Pc-scale jet propagation 

To our knowledge, first observational hints on jet ro- 
tation were provid ed by Echelle spec troscopy of knots 



in the HH212 jet (Davis et al. 2000). For one of the 



inner knots (SKI, 21)00 Au from the source), rotational 
velocities of 1.5km/s were derived, while other knots in 
the same jet/counter-jet did not follow the same sense of 
rotation (clockwise) in spite of the almost perfect side- 
symmetry of HH 212 and the fact the sense of r otation of 
the central molecular disk is clockwise as well ( Wiseman 
eraL][200ll). 



Davis et ai. ( 2000 ) estimate the angular momentum of 
an mner accretion disk of 500 AU. Assuming that all that 
angular momentum will be removed by a jet with mass 
ejection rate of 10% the accretion rate they conclude that 
HH212 may rotate with about 2km/s. If the disk angu- 
lar momentum is removed by viscous transport only, the 
jet rotational velocity would be lower. We find this ar- 
gument somewhat delusive as the high velocity jets is in 
fact launched close to the inner edge of the disk. Thus, 
the disk material which is ejected into the jet at these 
small radii has low specific angular momentum as it has 
lost most of its initial angular momentum during the ac- 
cretion process from larger to smaller radii. 

A correspon ding search for o utflow rotation in HH212 
in SiO shocks ( Lee et aL||2008 ) provided upper limits for 
rotational velocities only. Applying angular momentum 
conservation for steady-state the inferred jet launching 
radii can be constrained to < 1 AU. Most recent data of 
HH212 gave no clear sig nature of jet rotatio n, but also 
cannot exclude rotation ( Correia et al.J[2009 ). 

In the case of HH211, launching radii ot 0.06 — 0.15 AU 
were derived from possible jet rotation of 1.5km/s ( |Lee 
etal 1120071 [2009|. 



The NGC 1333 IRAS 4A2 system is another jet source 
for which recent SiO observations in dicate jet rotation 
by perpendicular velocity gradients (Choi et al. 2011). 
These features extend as far as almost 9000 A U from the 
central protostar with derived rotational speeds of about 
3km/s and in corotation with the underlying disk. 

2.2. Micro-jets from young stars 
The observed Do ppler-shift of emission lines in the 



DGTau micro-jet (Bacciotti et al. 2002) give sym- 



metric radial velocity ditterences of 6-15 km/s for the 
low-velocity jet component. It is difficult to imag- 
ine any other velocity pattern than rotation to explain 
these radial velocity gradient across the jet. The au- 
thors estimate the corresponding angular momentum 
flux in the low velocity component to J ~ 3.8 x 
lO-^Mmyr-^AUkms-^ 



More exampl es were detected by the same group ( Cof- 
fey et al.||20p 4 [), such as TH28 and RWAur or LkH^ 
l-br RW Aur Woitas et al. ( 2005 ) find toroidal velocities 
in the range ot 5 — 30kni/s at distances 20 and 30 AU 
away from the jet axis. Furthermore, both bipolar lobes 
rotate in the same direction. 

Another candidate source for jet rotation is HL Tau 
for which Fabry-Perot interferometry has indicated a 
radial velocity gr adient across the jet w ith a slope 
of 6km/s/arcsec (Movsessian et al. 2007). As for to 
DG Tau, it is the low velocity component which shows 
a significant radial velocity slope while the radial veloc- 
ity variation for the high velocity component is weak. 

2.3. Molecular outflows 

Indicate exists for rotation in molecular outfiows as 
wel l. From sub-mm da ta of the bipolar outflow Ori- 
S6 Zapata et al. (|2009) flnd bullets of in rotation of 
5 km/s, a CO shell around the jet rotating with 2 km/s 
(at 1000 AU from the outflow axis), and a CO envelope 
at distances of more than 2000 AU from the source still 
rotating with more then 0.5 km/s, and rotation observed 
up to 25,000 AU distance from the source. 

Another example i s CB 26, a molecular outfl ow from a 
young TTauri star (Launhardt et al. 2009). Detailed 



^ Observational indication of jet rotation had not been found by 
that time 



kinematic modeling accompanied by synthetic CO maps 
suggest rotational velocities of about Ikm/s 100 AU 
away from the outflow axis and at 1000 AU away from 
the source. Alternative scenarios such as a precessing 
jet entraining material are discussed for this source and 
seem feasible as well. 

2.4. Magnetic field indication in protostellar jets 

In general, there is indirect indication of magnetic fields 
being present in protostellar jets. One argument has 
always been the magnetic activity of protostars show- 
ing large-scale flares on radio and X-rays and features 
of dipolar accretion indicating on a kG stellar surface 
field strength and also on a dipolar field distribution (see 
(Andr e et al.][T988l [BOTlvier et al.|[2007| |Cabrit|[2007l )). 

Concerning the large-scale field structure, there is the 
unique example of TTau where cyclotron radiation have 
been detect ed which is aligned with the jet outfiow in 
this source (|Ray et al.|[T997). The infered field strength 



is in the range ot a G which is huge compared to what is 
required by current modeling. 
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However, recent observations of massive star forming 
regions have detected maser emission in jet sources which 
allows for a direct measurement of the magnetic field 
strength. One example is the detection of the 3D field 
structure aroun d the protostar CephAH W2 with mG 



field strengths (Vlemmings et al. 



emission indicating a (J. 2 mG held has been found re 



2010). Synchrotron 



lated to the knots in the jet source HH80-81 (Carrasco- 
Gonzalez et"aL||2010[ ). 



3. THE CASE OF EXTRAGALACTIC JETS 

Observational indication of rotation is lacking so far 
for extragalactic jets. In fact, we do not even have direct 
measurements of propagation speed for these sources. On 
the other hand, clear evidence exists for a helical mag- 
netic field structure in extragalactic jets indicated e.g. 
by gradients in t he rot a tion measure ac r oss th e jets (e.g. 
Gab uzda e t"al] (|2004P; iLyutikov et all (|2005P; ILaing et 



The existence of radio synchrotron emission along 
these jets on time scales beyond the synchrotron cool- 
ing time requires re-heating of the jet material, most 
probably by internal shocks. Together, the helical field 
structure and the internal shocks, provide the necessary 
prerequisites for the acceleration mechanism discussed 
below. Therefore, signatures of jet rotation can also ex- 
pected from these sources. 

In fact, there is one jet source where azimuthal jet mo- 
tion has been detected. Imaging interferometry of emis- 
sion lines have detected a braided structure in the jet 
system of NGC 4258, which could be niodeled by a tripl e 
helix of jet motion on 5-10 kpc scale (Cecil et al. 1992). 
This kpc jet is aligned with a nuclear VLBA radio jet 
(Herrnstein et al. 1997), thereby suggesting a physical 
connection. However, even if the helix model would ex- 
plain the physical conditions in the kpc jet, there seems 
to be no plausible explanation yet about the intrinsic 
process causing such a complex outflow structure. 

In this context it is interesting to note the relevance 
of jet rotation for modeling the heating of relativistic 
jet plasm a by shear and centrifugal acceleration for elec- 
trons (Rieger &: Mannheim|[2002 [Rieger fc Duffy||2004 | 
or protons ( [Dempsey fe Rieger||2009 ). 



4. MHD JETS AND ROTATION 

Magnetohydrodynamic jets are launched from rotating 
sources. The initial acceleration is by centrifugal forces 
along the magnetic field lines (magneto-centrifugal) . Be- 
yond the Alfven point along each field line, the outflow 
material becomes inertially so heavy that co-rotation of 
matter and magnetic field cannot longer be maintained. 
Thus, the field lines will be "bent" from a mainly poloidal 
into a helical geometry (a toroidal field component is in- 
duced). This allows the matter to "slide" along the field 
lines in toroidal direction. As a consequence, at large 
radii the jet material moves outwards in almost radial di- 
rection by conserving its total angular momentum. The 
angular velocity decreases with radius. 

Magnetohydrodynamic acceleration and coUimation of 
outfiows implies transfer of magnetic energy to kinetic 
energy, and, likewise, transfer of kinetic angular momen- 
tum in magnetic angular momentum. The exact balance 
between these quantities results from the local force- 
balance throughout the outflow. The solution of the 



equations is complex and requires semi-analytical mod- 
eling and/or numerical MHD simulations. 

For understanding the considerations of this paper, it is 
helpful to re-call some basic properties of axisymmetric, 
stationary-state theory of MHD winds and jets (see e.g. 
,Pudritz et al.,(2007) ). In general, steady-state MHD im- 
plies the existence of conserved quantities along the mag- 
netic flux surfaces. In the case of axisymmetry flve such 
quantities can be deflned. The flrst one is the magnetic 
flux, ^(r, z) = (l/27r) J BpdA, indicating that magnetic 
field lines lie on magnetic flux surfaces. The second one 
is the mass flux 77(5") = / pVpdA per magnetic flux sur- 
face ^P. The third one is Ferraro's iso-rotation parameter 
i7F(^) = {v-pB^/B-p — v^)/r, often interpreted as angu- 
lar velocity of the field lines. The fourth one is the total 
(specific) angular momentum 



L(*) = rv^ 



rB^Bp 
AirpVp 



ilpr%. 



The fifth one is the total (specific) energy 



1 2 

2"P 



rB^Bpflp 
AnpVp ' 



(1) 



(2) 



where, for simplicity, we have omitted gravity and gas 
pressure (cold wind assumption). In steady-state theory 
E{^>), L{'^), and 77 (^') are governed by the regularity 
condition across the fast magnetosonic point, the Alfven 
point at r — rj^, and the slow magnetosonic point, re- 
spectively. 

Essentially, both the total energy and the total angu- 
lar momentum (flux) consist of a magnetic and a kinetic 
contribution. A classical Blandford-Payne-type outflow 
is launched with » Vp and i?0 < -Bp. Beyond the 
Alfven point decreases and the magnetic held becomes 
dominated by B^ > Bp. Likewise, kinetic angular mo- 
mentum of the ejected jet material, which is co-rotating 
with the disk, is increasingly transferred to the outflow 
magnetic angular momentum. 

The existence of conserved quantities along the field 
lines can be used to physically link the very inner region 
where the jet is launched (not yet possible to resolve 
observationally) to the asym ptotic region of the outflow 
(accessible for observations). [Anderson et al.| ([2003) have 
applied this formalism to the observed rota tional fea- 
tures in the DGTau jet ( Bacciotti et al.|2002[ ), obtaining 
launching radii of 0.3-4.0 AU. While this approach could 
be applied to the low velocity jet component, it does not 
fit, however, the high-velocity jet component of DGTau. 
Using a simil ar method for the bipolar jet from RW Aur, 
Woitas et al.|([2005) find launching radii less than 0.5 AU 
and a magnetic lever arm (Alfven radius to launching 
radius) of about four. 

In the following we will argue that the classical station- 
ary MHD wind approach discussed above is questionable 
as sole explanation of rotating outflows. In the next sec- 
tion we will suggest and discuss an alternative scenario. 
Our main argument is that radial velocity measurements 
of stellar jets and outflows usually consider line emis- 
sion of shocked gas or molecules. For example, [Lavalley^ 
Fouguet et al. ( 2000 ) prove shock excitation of f orbid- 
den emission lines m DG Tau and Woitas et al. ( 2005 1 
explicitely state that rotation is more evident in the tor- 
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bidden lines [01] and [Nil]. During shock transition, the 
jet material changes its dynamical state instantly. Shock 
transitions are not included with the classical stationary 
state MHD wind theory. Instead, conservation laws must 
be applied in form of jump conditions, and the direct con- 
nection of the observed rotating material with the disk 
by ideal MHD conservation laws is questionable. This ar- 
gument may hold in particular for the large scale jets and 
outflows as e.g. for IIII212 or molecular outflows. The 
molecular material might also trace entrained interstellar 
material. In this case, it has not been launched from the 
central region with intrinsic angular momentum. 

Early work on steady-state MHD jets has suggested 
that jets do expand substantially from their launch- 
ing area close to the star to the asymptotic (observed) 
regime. Expansion factors of 100-1000 (from 0.1 AU up 
to 100 AU jet diameter) have been suggested. These clas- 
sical solutions predict Alfven radii i?A of about ten foot 
point radii, Ra — lORp and, therefore, a substantial de- 
crease of the toroidal outflow velocities for radii r » Ra 
to v^/vjet < 0.001. This is due to the large initial open- 
ing angles whi ch are a resu lt from non-self-similar MHD 
solutions (see Fendt et al. (1995); Fendt & Camenzind 
([1996), but als o recent smiuiations ot ste llar wind mag- 



(2008 2010) resuhing in a 



netospheres by Matt et al. 
similar field geometry). 

We conclude that while being consistent with station- 
ary state MHD jet formation, the formalism suggested by 
Andersen et al. and subsequently applied by other au- 
thors is probably not a fully realistic picture to explain 
the observed jet rotation. 

5. MHD SHOCKS CRANKING JETS 

In the following we discuss another framework of how 
jets and outflows may be set in rotation. The basic idea 
is that a compression of the toroidal magnetic field com- 
ponent in MHD shocks results in a Lorentz torque in 
toroidal direction which accelerates the jet material. The 
jet kinetic angular momentum is retrieved from the mag- 
netic angular momentum content of the jet by the shock. 
This mechanism works only for helical magnetic fields. 
While in stationary state MHD the kinetic and magnetic 
angular momentum is gradually exchanged along the flux 
surfaces, the shock wave leads to a sudden exchange. 

5.1. Lorentz torque 

The MHD Lorentz force is Fl - J x B - ( V x B) x B. 
Jet rotation can be achieved by its toroidal component 

(3) 

Across the shock the toroidal field component is com- 
pressed, thus giving rise to a toroidal force component. 
The compression ratio itself depends on the jet dynamics. 
For simple geometries the field compression can be esti- 
mated by the Rankine-Hugoniot jump conditions. More 
complex configurations require numerical simulations to 
follow the shock evolution. 

The magnetic fiux is conserved across the shocl|^ 
[Bz]i = 0, but the perpendicular magnetic field compo- 
nent typically undergoes a change. In case of cylindrical 

^ Index '1' denotes the upstream region and index '2' tlie down- 
stream region 



- (V X B^) X Bp. 



coordinates (which we are using here) , the iJ^-component 
replaces the perpendicular component in Cartesian coor- 
dinates, resulting in an expression for the compression 



Bi 



2mt 



■/Bl 



r *i2 
[P]i 



Bl 



(4) 



( Uchida et al.| |1992), where = pv^ denotes the (con- 
served) axial rnomentum, (p*) = 0.5 (l/pi-|-l/p2) the av- 
eraged inverse density from upstream and downstream, 
[p*]^ the difference in the inverse density upstream and 

downstream, and {B^)^ = {B)^ — B^ the average toroidal 
field strength. 

The dynamical structure of the propagating shock 
could be complex. As typically for this setup, four shocks 
are existent, accompanied by a contact discontinuity, 
each sa tisfying the jump cond itions (see e.g. l Uchida et al.| 
(| 1992 1 ; |Ryu fc Jones] |l995[ )). Considering the toroidal 
field strength, an increase within the first MHD (fast) 
shock is followed by a decrease across the second MHD 
(slow) shock, then again by an increase in the third (slow 
MHD) shock, and another decrease in the fourth (fast 
MHD) shock ( going from dow n stream to ups tream; see 
below and , Uchida et al.| ^1992[ ); |Ryu fc Jones] ( |1995J ). 

A vast number ot multi-dimensional jet propagation 
simulations have been published. Only few of them apply 
a helical magnetic field configuration and some of those 
have already provided numerical indication for toroidal 
velocity changes due to toroidal magnetic torques in jet 
shocks. A connection to or predictions of the observable 
jet rotation has not been made (such observations were 
not available in that time) . Ear ly simulations of jet prop- 
agation do not show rotation (Clarke et al. 1986) since 
treating a purely toroidal magnetic held only! To our 



knowledge, the MHD simulations by Kossl et al. ( 1990 1 



were first indicating rotational velocity changes caused 
by a toroidal magnetic torque across a shock. Applying 
jet parameters such as Mach numbers of 4, plasma- /3 = 1, 
and a helical field configuration B, 
velocities obtained were rather small vs < 0.1c, 



Bp the rotational 



Uchida et al] ( |1992[ ); |Todo et al | ( [1992) present MHD 
simulations in 1.5D and 2.5D which m particular com- 
pare 1.5D simulations in cylindrical coordinates to ana- 
lytical solutions of steady MHD shocks. In their case, the 
rotational velocities obtained for Ma — 1.0, 2.5, 6.0 simu- 
lations are comparatively high - in fact up to 10-40speed, 
and, thus, not realistic compared to the observed values. 
In fac t, in a forthcoming paper on 3D simulations |Todo] 



et al. (1993) discuss such high azimuthal velocities as 
fo remost indicator of a jet instability. 

( 1999 ) investigate jet stability under the 
ical magnetic fields, but do not discuss jet 
rotatio n in particular. Simulations by IStone fc Hardee] 



Rosen et al. 
influence ot he 



(20001 consider both steady and pulsed MHD jets with 
helical magnetic field, in particular their cooling and sta- 
bility properties, but again do not stress jet rotation. 

5.2. Angular momentum budget 

A basic prerequisite for the applicability of the process 
discussed in this paper is a sufficient amount of magnetic 
angular momentum available, which can be converted 
within the shock region. We may estimate the angu- 
lar momentum budget from numerical models of the jet 
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formation region which provide solutions to the MHD 
equations along the accelerating and expanding jet. 

The amount of magnetic angular momentum flux Lmag 
in respect to the kinetic angular momentum flux Lkin de- 
pends crucially on the MHD stat e of the flow. Analytical 
estimates (Blandford fc Payne 1982[ IPelletier & Pudritz 

c n( 



Table 1 

Typical astrophysical scaling parameters for jets from young stars 
and AGN 



ows L 



1992 1 show that for sub-Alfvemc 
wtiile for marginally super-fast flows Lmag 



mag 



kin; 

2Lkin- For a 

(poloidal) fast magnetosonic Mach number Mfm = 2 the 
kinetic angular momentum flux dominates, imag/ikin — 
2/M|m=0.5.^ 

Numerical simulations of MHD jet formation from Ke- 
plerian disks show, however, that typically only the in- 
ner jet (the pa rt close to the jet axis) becomes super-fast 
magnetosonic ( Ouyed fc Pudritz 1997; Krasnopol sky et 

«i innn, t7,,^ jj^^^^^T^^k:^ or^fio. w...^A^ . ij^^^C^^^ 



R. 1999; Fe ndt fc C'emel jic 2002; Fcndt 2006; PortlTfc 
^endt 2010|. This suggests an asymptotic jet structure 
consisting of a narrow core dominated by kinetic angu- 
lar momentum, which is wrapped by an envelope dom- 
inated by magnetic angular momentum flux. It is this 
envelope around the core jet, which i) (shock-) interacts 
with the ambient medium, and ii) is accessible to the 
observations. Thus, for the purpose of this paper it is 
convenient to assume that a substantial (i.e. sufficient) 
amount of magnetic angular momentum is present in the 
asymptotic jet. 

Numerical simulations of MHD jet formation from Ke- 
plerian disks also show quite a variety in the asymptotic 
jet outflow they produce (and which is su bject to the 
present paper). We refer to our past work ( Fendt|[2006 l 
which in very detail investigated how the launching con- 
ditions as the disk magnetic held profile and strength and 
also the mass flux distribution affects the asymtotic flow. 
In particular the radial profile of the disk wind magne- 
tization (steep or flat) determines the collimation of the 
outflow (weakly, resp. strongly coUimated). We there- 
fore decided to scan a rather large parameter range in 
the typical jet characteristics as Mach number or plasma 
beta. 



Recently Ramsey fc Clarke (20111 have performed 
jet formation simulations spanning from the disk surface 
well into the asymptotic regime, enclosing a huge, yet 
unprecedented area of 256 x 4096 AU size. These simula- 
tions show i) that a well ordered, global held distribution 
remains present also for the most distant outflow regime, 
ii) that the asymptotic flows launched with different ini- 
tial plasma-/? tend to approach asymptotic states with 
average plasma-/? close to unity (although there is still 
quite a spread during the simulation time scales consid- 
ered), and iii) that flows of different initial plasma-/? may 
reach quite different asymptotic states concerning mass 
fluxes of velocities. 

In summery, we decided to run our simulations for a va- 
riety of leading jet parameters, considering also the fact 
that jets from protostars and AGN could be intrinsically 
different. 

6. MODEL SETUP FOR THE MHD SIMULATIONS 

In order to quantify the toroidal torque in jets and the 
subsequent toroidal velocity structure, we perform MHD 
simulations of jet propagation applying different flow ge- 
ometries, magnetic field distributions, and jet dynamical 
parameters. In order to demonstrate the efficiency of 
the proposed mechanism, we inject the jet material non- 





Stellar jets 


AGN jets 


Viet 


~ 300km/s 


~ 0.3c 


Pjct 


~ 10*cm-3 


~ 10-*cm-3 




10 - lOO^G 


10-100 G 



rotating, thus with vanishin g kinetic angular mo mentum. 

We use the PLUTO code |Mignone et al.|2007| ) to solve 
the standard set of ideal MHD equations. We apply an 
ideal equation of state with polytropic index 7 = 5/3. We 
also discuss one example simulation for which cooling is 
considered. 

In order to re-scale normalized code variables to astro- 
physical magnitudes, one may apply typical numbers val- 
ues for protostellar jets or AGN jets as shown in Tab. [ij 

6.1. Simulations in 1.5D 

With "1.5D" we denote axisymmetric simulations tak- 
ing into account all three spatial coordinates of the vector 
fields, however jPutting one of the vector field component 
to zero initialljp We also limit the grid extension to only 
one cell in this airection. 

The general setup is that of a typical MHD shock tube 
test simulation. High-density/high pressure gas of high 
velocity is continuously injected from the left-hand side 
of the box (considered as the jet flow). Low-density/low 
pressure material is placed at the right-hand side of the 
box (considered as the ambient gas) and is allowed to 
move out of the box. With this setup we extend pre- 



viously exis ting simulations ( [Uchida et al.||1992| Ryu & 
Jones 1995 ) to the specific question of jet rotation and 
for parameter studies of that. 

We apply different magnetic field configurations. For 
the simulations in cylindrical coordinates we assume an 
initial axial field Bz{z) which is constant along the grid. 
The initial toroidal field B^{z) is either only present on 
the left side (upstream) or constant along the grid. For 
comparison with the literature and in order to test the 
1.5D approach in cylindrical coordinates, we have com- 
puted the same problem also in Cartesian coordinates 
and have found perfect agreement (see appendix). 

The numerical resolution varies from 5000 - 20000 grid 
cells for 20 physical scale lengths. 

6.2. Simulations in 2.5D 

With "2.5D" we denote axisymmetric simulations tak- 
ing into account all three vector components. By def- 
inition the simulations sustain axisymmetry as the </>- 
derivatives are not considered in this option of the code. 
A jet of unity radius is injected into a box of ambient gas 
with constant density, pressure, and longitudinal mag- 
netic field. We apply a non-equidistant grid with high 
resolution across the jet beam and the interaction region 

* That is the r-component in cylindrical coordinates. We have 
also performed comparison simulations in Gartesian coordinates 
with an initially vanishing x-component. For both cases we have 
approved that these components remain negligible during the sim- 
ulation, i.e. < 10~^*. 
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with the ambient gas (r < 1.5, 100 cells), and a some- 
what lower resolution in the ambient gas (1.5 < r < 7.0, 
200 cehs). 

It would be interesting to extend the present investiga- 
tion to the 3D case in a future paper. Although the prin- 
ciple mechanism proposed here will be the same in the 
3D case, the overall jet structure could be different and 
also the characteristic instabilities affecting the jet flow. 
In particular large-scale instabilities as the m = 1 kink 
mode could appear and disturb the jet flow. Such insta- 
bilities would develop beyond the Alfven surface (thus in 
the super- Alfenic flow), and may de stroy the large-scale 

This is a serious 



2009) 



helical field structure ( Moll 

argument, ho wever, it aftects the st andard model of jet 
rotation (e.g. Anderson et al. (20031) even more. On the 



other hand 3L) simulations of supersonic have revealed a 
quite similar general structure compared to the 2D case 
( [O'Neill et "aL||2005|). Also, 3D simulations of jet for- 
mation by Ouyed et al. ( 2003 ) revealed a self-stabilizing 
ability of the jet due to its axial spine. We finally note 
the fact that we do observe large-scale collimated out- 
flows which tells us that jet stability must be maintained 
by some (yet unknown) process. 

6.3. Main simulation parameters 

Our simulations are characterized by the following set 
of parameters. 

We apply different jet densities pjet, while the initial 
ambient density is always the same, pext = 1.0. The jet 
inflow velocity is ujet- In the case of 2.5D simulations 
this is the maximum jet velocity as we prescribe a cosine 
proflle for the injection speed across the jet inlet. The 
ambient gas is initially at rest. The jet gas pressure is 
Pjet and the ambient gas pressure is Pext- The strength 
of the magnetic field is parameterized by the (squared) 
jet internal poloidal Alfven Mach number. 



M 



Ap.jct 



and the (squared) jet internal total Alfven Mach number. 



Ma. 



jet 



■°P,jct 



With the (squared) jet internal sound speed jet — 

7-fjct/Pjet and Mach number Ms jot = Wpjot/csjot the jet 
internal fast magnetosonic Mach number is 



Mpjet = I M 



AJct 



SJct 



-1/2 



The external Mach number relates the poloidal jet veloc- 
ity to the sound speed in the external medium, Mg oxt = 

'^pjet / Cs,cxt • 

Note that the Mach numbers are defined in respect 
to the poloidal velocity, considering MHD waves in the 
jet propagation direction. Tables [2] [4] summarizes the 
parameter setup for the 1.5D simulations, whereas Tab. [3] 
shows the parameter space for the 2.5D simulations. 

The initial magnetic field _Bz-component is constant 
over the whole computational domain. We run simula- 
tions for which the toroidal field component is in- 
jected into the ambient gas by the jet and others where 



a shock propagates along an already existing helical field 
structure. For the jet toroidal field we have applied a lin- 
early increasing radial profile and a sine-profile. While 
the linear profile is force-free within the jet, we find the 
the sharp cut-off at the jet boundary leads to current 
sheets which disturb the flow ab initio and could not be 
accepted. We therefore discuss the simulations applying 
the (non force-free) sine profile only. 

7. RESULTS AND DISCUSSION 

We first discuss simulations resulting from a 1.5D 
setup. These can be considered as toy simulations of ax- 
isymmetric shocks in a helical field and test cases of the 
model. We then continue with 2.5D axisymmetric sim- 
ulations of jet propagation with parameter ranges more 
comparable to observational values. 

7.1. Jet shock propagation in 1.5D 

Quite some literature exists on one-dimensional MHD 
shock simulations, as this setup provides one of the stan- 
dard tests for MHD codes. In particu lar we refer to 



'Uchi da et aL] ( jl992[ ); 
repeat and extend a 



Ryu fc JonesMl995[ ). We decided to 
:ew of these 1.5D simulations for the 



following reasons. Firstly, they serve us as test case for 
our code for the 1.5D setup in cylindrical coordinates as 
we can directly compare them to simulations in the lit- 
erature in Cartesian coordinates which are widely used 
also for analytic considerations. Secondly, this simple 
setup will most clearly demonstrate the essential mecha- 
nism leading to a toroidal jet acceleration across a shock. 
Thirdly, it allows us easily to explore the typical pa- 
rameter space for protostellar jets (see Tab. [2|. In the 
majority of our simulations the jet is injected as non- 
rotating. This simplification is made in order to demon- 
strate the acceleration mechanism by the shock, and is 
in general only valid for a certain iso-rotation parameter 
= (vM/iBpr). 

We have performed several test simulations, of which 
two are shown in the appendix. That is simulation 
RJll repeating the simulation by jRyu & Jonesj (|1995| 
(see their Fig. 2a), done as well in Cartesian coordinates 
(Fig. 14 1 . Simulation RJ09 is an extension of this param- 
eter setup using however cylindrical coordinates (Fig. 15 1 
and considering only one perpendicular vector compo- 
nent (no radial magnetic field, no radial motion). 

Our first science example is a shock moving along a 
he lical magnetic field (run RJ17), similar to the setup 
of |Ryu fc Jones (1995). This is comparable to a jet in- 
ternai shock propagating within the jet channel, caused 
e.g. by a sudden increase of jet injection velocity. Figure 
[l] shows the dynamical state after two dynamical time 
scales. The typical shock structure is visible, consisting 
of four shocks which propagate back resp. forth in re- 
spect to the contact discontinuity. The hydrodynamical 
state of the different shock regions varies, each satisfying 
the Rankine-Hugoniot jump conditions. The toroidal ve- 
locity panel in Fig. [T] clearly shows how the material is 
set into rotation when crossing the shock fronts. 

Across the shock the toroidal field becomes com- 
pressed, and the steeper field gradient implies a toroidal 
torque on the shocked material. The material down- 
stream and upstream the shock rotates in opposite di- 
rection, seen also in axial profile of kinetic angular mo- 
mentum. This is a consequence of a change of sign 
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Table 2 

Parameter of axisymmetric 1.5D simulations. The external density is pcxt = 1-0 for all runs. The toroidal velocity is given before and 
after the shock front v^ -i,v^ 2- Remarks: UT05 follows the same parameter setup as UT04, but c onsid ers cooling. UT06 follows the same 
parameter setup as UT04, except for the negative rotation velocity. For comparison with Figs. [T|2] note the PLUTO magnetic field 

normalization B —> B/\/4tt. 



ID 


Pjot 




"■/..jet 




Pcxt 




-B*,ict 


.B^.cxt 


A/Ap,jct 


AfA.Jct 


A'^Fjct 


A/s.cxt 




RJ03 


2.5 


1.0 


0.0 


1.53 


0.27 


0.11 


1.70 


0.0 


50.0 


3.5 


0.88 


1.5 


0.0/-0.163 


RJ04 


2.4 


1.0 


0.0 


0.32 


0.27 


0.11 


1.56 


0.0 


50.0 


3.5 


1.82 


1.5 


0.005/-0.235 


RJ06 


2.4 


1.0 


0.0 


0.026 


0.27 


0.18 


1.00 


0.0 


30.0 


5.5 


4.42 


1.5 


0.017/-0.22 


RJ17 


1.08 


1.2 


0.0 


0.95 


1.0 


2.0 


2.0 


2.0 


2.21 


1.56 


0.84 


0.93 


+0.091/-0.088 


RJ18 


1.08 


1.2 


0.0 


0.95 


1.0 


2.0 


2.0 


0.0 


2.21 


1.56 


0.84 


0.93 


+0.084/-0.31 


RJ19 


1.28 


1.5 


0.0 


1.15 


1.0 


2.0 


2.0 


0.0 


3.01 


2.13 


1.06 


1.16 


+0.078/-0.29 


RJ21 


1.28 


1.2 


-0.2 


0.95 


1.0 


2.0 


2.0 


0.0 


2.41 


1.70 


1.10 


0.93 


-0.12/ -0.41 


RJ22 


1.28 


1.2 


0.0 


0.95 


1.0 


0.5 


2.0 


0.0 


9.63 


2.34 


0.98 


0.93 


+0.019/-0.33 


RJ23 


1.28 


1.2 


0.0 


0.95 


1.0 


2.0 


0.5 


0.0 


2.41 


2.34 


0.98 


0.93 


+0.029/-0.078 


UTOl 


1.0 


4.0 


0.0 


0.1 


0.01 


4.0 


2.4 


2.4 


3.55 


3.04 


2.90 


31.0 


-I-0.7/-0.7 


UT02 


1.0 


4.0 


0.0 


0.1 


0.01 


4.0 


2.4 


0.0 


3.55 


3.04 


2.90 


31.0 


-I-0.7/-0.4 


UT03 


1.0 


4.0 


0.0 


0.1 


0.01 


1.0 


2.4 


0.0 


14.2 


5.44 


4.76 


31.0 


-I-0.14/-0.54 


UT04 


1.0 


4.0 


0.1 


0.1 


0.01 


1.0 


2.4 


0.0 


14.2 


5.44 


4.76 


31.0 


-t-0.24/-0.48 


UT05 


1.0 


4.0 


0.1 


0.1 


0.01 


1.0 


2.4 


0.0 


14.2 


5.44 


4.76 


31.0 


4-0.25/-0.46 


UT06 


1.0 


4.0 


-0.1 


0.1 


0.01 


1.0 


2.4 


0.0 


14.2 


5.44 


4.76 


31.0 


+0.043/-0.58 



in the toroidal Lorentz force ^ BzdBtp/ dz across 
the shock. With increasing time, the shock evolution 
demonstrates the continuous conversion of magnetic an- 
gular momentum flux into kineti c angu lar momentum 
flux (not shown here, but see Figs. 12|13 below). The the 
domain-integrated kinetic angular momentum increases 
in time. Due to the downstream boundary condition for 
the toroidal magnetic field, the net amount of kinetic an- 
gular momentum along the flow remains zero: Implicitly, 
the outflow conditions for angular momentum is kept the 
same as the inflow condition. 

Our next example RJ18 applies similar flow parame- 
ters, but corresponds to a model setup where a jet is 
injected into an ambient medium threaded by an axial 
magnetic field only (the ambient toroidal magnetic field 
is set to zero). Figure [2] shows the dynamical state at the 
same time as for RJlT. The difference comes from the 
magnetic angular momentum flux which now decreases 
along the flow, while being transferred to kinetic angular 
momentum. For the parameters applied the toroidal ve- 
locities obtained are about 25% of the jet velocity, and 
much higher than the observed values. 

Simulations RJ17 and RJ18 exhibit a similar hydrody- 
namic shock structure (considering wjot, p, P)- Because of 
the stronger toroidal field coiling, the MHD torque and, 
thus, the resulting toroidal velocities are increases by a 
factor of three. 

Figure [3] shows simulation RJ21 in which a jet rotating 
with w^jot = —0.2 is injected into the ambient medium. 
By choice, the injected toroidal field direction is anti- 
aligned with the toroidal velocity. The shock structure 
of RJ21 is very similar to RJ18 which has otherwise the 
same hydrodynamic and magnetic parameters. Compar- 
ing RJ18 to RJ21, the resulting maximum jet rotational 
velocity is enhanced from — —0.31 to v^p — —0.41, 
while the rotational velocity jump across the receding 
shock is the same in both cases, Aw^ = 0.74. — —0.41, 
but also from = —0.31 to = —0.41. Figure [3] 
also shows the hydrodynamic and the magnetic angular 
momentum contribution to the total angular momentum 



budget (see Eq. [T] but note that we measure the a.m. at 
unity radius r = 1). For each of the four shocks the to- 
tal specific angular momentum is conserved in the shock 
frame satisfying the Rankine-Hugoniot conditions. It can 
be clearly seen that the magnetic angular momentum de- 
creases across the shock as the kinetic angular momen- 
tum contribution increases. 

It is also interesting to compare simulations RJ22 and 
RJ23 (see Figs.[4j[5]). Both start off with the same initial 
magnetic angular momentum budget and Foynting flux 
^ BpB^, however, RJ22 has B^, — 4i3p, while for RJ23 
it is Bp — ABfj,. From this, one may expect also a simi- 
lar toroidal magnetic torque and, thus, a similar toroidal 
velocity, in particular since a similar compression rate 
in the hydrodynamic parameters and also in B^ applies. 
However, we find that the resulting maximum toroidal 
velocities substantially change from 0.19 to 0.29 for the 
receding shock, or from -0.33 to -0.078 for the preceding 
shock and similarly for the kinetic angular momentum 
(see Fig. B. The reason is that the stronger poloidal 
field in RJ22 broadens the shock structure and thus dis- 
tributes the kinetic angular momentum over a wider spa- 
tial area. This is a typical property of C-shocks. Note, 
however, that the total, volume-integrated kinetic angu- 
lar momentum in the shocked region is the same since 
the volume of the shocked gas differs (the shock widens 
from Ar = 1.5 to 3.5 for the preceding shock). 

We have also performed a comparison simulation with 
applying a simple power-law cooling function. Our simu- 
lations indicate a minor effect of cooling on jet rotation. 
The shock width decreases slightly, as the gas pressure 
decreases by cooling (this is a 10%-effect for the chosen 
cooling function) . The gas density and toroidal magnetic 
field strength increase resulting in a slightly increased 
maximum toroidal velocity. These preliminary results 
need to be further investigated in a future paper. 

In summary, our 1.5D simulations clearly show how a 
jet flow in a helical magnetic field is forced into rota- 
tion, when the toroidal fleld component is compressed 
across the shock, thus exerting a Lorentz force on the 
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Figure 1. Axisymmetric 1.5D shock simulation RJ17 in cylindrical coordinates. Shock evolution at dynamical time t = 2 with a grid 
resolution of 1.000 cells / unit length. Shown is density; gas pressure; axial, and toroidal velocities; axial, and toroidal magnetic field 
strengths, and the kinetic angular momentum (at r = 1) (from top left to bottom right). 



jet material. Typically, jet rotation velocities are about 
3 — 10% of the jet axial velocity for Alfven Mach numbers 
MAp,jet :^2-3 (see Tab. [2|. 

7.2. Jet shock propagation in 2.5D 

We now discuss axisymmetric simulations taking into 
account all 3 vector components. This approach allows 
to investigate also the lateral expansion of the jet and 
its interaction with the ambient medium. We have run 
a parameter study, mainly considering different values 
for the jet magnetosonic Mach numbers and plasma beta 
(see Tab. [3]). 

The general boundary condition applied for the 
jet toroidal magnetic field follows a profile — 
S0,max(sin((2.Or - 0.5))7r -I- 1.0) for r < 1.0, and van- 
ishes for larger radii. This is not a force- free con- 
figuration, however, it avoids the strong current sheet 
usually launched by sawtooth-like profiles. The initial 
poloidal field distribution is constant and purely longi- 
tudinal. The radial profile of the jet injection velocity 
follows a cosine profile, Vinjir) = Vjet cos(r7r/2) for r < 1. 
We have also run simulations with a time-dependent 
injection velocity, Vjot{t) — 1.0 + 0.5 * sin(0.3i), thus 
Vm3it,r) — Winj(r)*(1.0-|-0.5*sin(0.3i)), or, alternatively 
with a longer period ~ (1.0 + 0.5 * sin(O.li)). 

Simulation run ROl considers a high external Mach 
number Mg^oxt = 10, a comparatively weak jet magnetic 
field corresponding to high magnetosonic Mach numbers, 
and a high plasma beta - thus, a jet rather close to the 
hydrodynamical limit. Figure [6] shows the time evolution 



of the toroidal velocity of this jet for the dynamical time 
steps t = 10, 30, 40, and the corresponding toroidal mag- 
netic field strength, axial velocity, and toroidal torque for 
time t = 30. This images show the inner domain r < 5.0 
of the simulation box r < 7.0. 

The simulation shows the rotating material surround- 
ing the inner jet beam in a cocoon-like structure. 
The injected jet material is non-rotating (by defini- 
tion of the boundary condition). For run ROl the ob- 
tained rotational velocities are comparatively low v^f, < 
10^'^Ujct, resulting from the low magnetic field strength 
in spite of the high shock compression ratio in the near- 
hydrodynamical limit. 

During the initial time steps the termination shock is 
rotating with positive sense, the bow shock with negative 
sense. The toroidal magnetic field compression is clearly 
seen in the termination shock, but also in surrounding 
cocoon. The cocoon toroidal field is enhanced partly 
by the radial expansion of the shocked gas, partly by 
winding up the poloidal field by the rotating jet material. 

The map of the toroidal torque demonstrates the in- 
terrelation between Lorentz force and rotation. The jet 
head is rotating in positive direction {t = 40, grey black- 
ish colors) with the toroidal torque in the same direction 
(white colors). The rotation velocity measured at a cer- 
tain time is, however, resulting from the time-integrated 
torque on a parcel of material along its path. 

Simulation run R12 considers a high (external) Mach 
number jet with weak magnetic field, similar to ROl, but 
representing a less extreme case (see Fig.[7]). The toroidal 
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Figure 2. Axisymmctric 1.5D simulation RJ18 in cylindrical coordinates. Shock evolution at dynamical time t = 2 with a grid resolution 
of 1.000 cells / unit length. Shown is density; gas pressure; axial, and toroidal velocities; axial, and toroidal magnetic field strengths, and 
the kinetic angular momentum (at r = 1) (from top left to bottom right). The toroidal magnetic field in the ambient medium is vanishing 
(compare to simulation RJ17 in Fig.jlJ. 



velocities induced range from -0.07 to 0.033 correspond- 
ing to 1-3 %. This is in the range of the observationaUy 
indicated velocities. Also the jet dynamical parameters, 



M 



A,jct 



30, Ms, 



5, are observationaUy indicated. 



Simulation run R13 considers the same parameters ex- 
cept a lower plasma-/3, thus a stronger poloidal magnetic 
field. 

For jets with lower Mach numbers (e.g. R09, RIO) even 
higher toroidal velocities can be produced. The main 
contribution to the rotating material is from entrained 
(and shocked) ambient gas. 

Simulation R16 runs with the same parameter setup 
as R12, but in difference applies a time-dependent injec- 
tion speed, Vjct{t) — Wjct,o(l-0 -I- 0.5sin(0.3t)). This leads 
to subsequent internal shocks along the jet, which due 
to B^-compression also force the jet material in rotation 
(see Fig. [8j|9]). In addition, each of the subsequent shock 
fronts further entrains ambient material. Accompany- 
ing shocks set this material into rotation as well. For 
this setup, it might be interesting to estimate the astro- 
physical parameters. Assuming a jet radius of 20 AU, 
the separation between the generated knots would corre- 
spond to about 600 AU (one grid length) . The knot pat- 
tern separation would then correspond to a time scale 
of 15 years, estimated from the 50 time units the knot 
pattern needs to traverse the 600 AU grid (one time unit 
At — 1 corresponds to about lO^s assuming a jet mean 
velocity of 300km/s). This time scale is of course put in 
by hand, but is not very different from observationaUy 



derived numbers. 

Simulation run R19 with lower Mach / Alfven Mach 
numbers compared to R16 has a slightly higher frequency 
in the injection speed variation. Together this results 
in more internal substructure , se en e.g. in the toroidal 
velocity distribution (see Fig. 10 1, but also in the other 
dynamical variables (not shown) . (and density, |6] not 
shown) . 

7.3. Oblique reflection shocks 

Jets with low (internal) Mach numbers develop stand- 
ing oblique relection shocks internal to the jet flow which 
do not further move as jet and bow shock propagate. As 
an example Fig. [TT] shows the dynamical state of sim- 
ulation run R12 at an intermediate time step t — 50. 
The reflection shocks are most clearly seen in the den- 
sity distribution and the radial velocity distribution. The 
change in propagation speed along the jet axis is about 
15% across each shock region. 

Along the jet flow the toroidal Lorentz force is predom- 
inantly negative leading to an according rotation of the 
jet material. The shock structure consists of a thin layer 
with positively directed toroidal Lorentz force followed 
(direction downstream) by a extended region of posi- 
tively directed toroidal Lorentz force. The latter force 
is dominating leading to an overall negative rotation of 
the jet beam. 

Comparing the toroidal Lorentz force and toroidal ve- 
locity distribution in Fig. 11 gives interesting insight on 
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Table 3 

Parameter study of axisymmetric 2.5D jet simulations. The ambient density is poxt = 1-0 for all simulations. The minimum and 
maximum values for the toroidal velocity < f ^ jot,(min,max) > ^''s derived from a histogram statistics, i.e. indicating the maxima for the 
bulk toroidal motion (there is always a tail in the toroidal velocity distribution of even higher velocities). The toroidal velocity maxima 
are estimates when the jet flow has penetrated most of the computational domain, but its bow shock has not yet left the grid. 



ID 


Pjot 




-Pjct 


-Poxt 


MApJot 






/3 


< 'y,^,jct,(mi„,max) > 


R,cmarks 


ROl 


5.5 


1.0 


0.006 


0.0015 


2000.0 


300.0 


10.0 


20.1 


-0.00014. 0.000065 




R02 


0.5 


1.0 


0.002 


0.006 


200.0 


30.0 


10.0 


10.1 


-0.0023, 0.04 




R03 


1.0 


1.0 


0.015 


0.006 


200.0 


30.0 


10.0 


10.1 


-0.028, 0.023 




R04 


1.0 


1.0 


0.012 


0.024 


200.0 


30.0 


10.0 


5.0 


-0.022, 0.025 




R05 


1.0 


1.0 


0.9 


0.15 


20.0 


3.0 


10.0 


2.1 


-0.42, 0.25 




R06 


1.0 


1.0 


0.9 


0.15 


50.0 


3.0 


10.0 


1.0 


-0.67, 0.40 




R07 


1.0 


1.0 


0.9 


0.15 


50.0 


3.0 


10.0 


2.0 


-0.55, 0.40 




R08 


1.0 


1.0 


0.5 


0.15 


50.0 


3.0 


0.5 


2.0 


-0.22, 0.15 




R09 


1.0 


1.0 


0.5 


0.15 


50.0 


3.0 


0.5 


2.0 


-0.15, 0.07 




Rll 


1.0 


1.0 


0.1 


0.15 


20.0 


3.0 


2.0 


1.1 


-0.18, 0.05 




R12 


1.0 


1.0 


0.001 


0.024 


200.0 


30.0 


5.0 


1.1 


-0.0023, 0.015 




R13 


1.0 


1.0 


0.0001 


0.024 


200.0 


30.0 


5.0 


0.1 


-0.0031, 0.010 




R15 


1.0 


1.0 


0.01 


0.067 


50.0 


10.0 


3.0 


1.1 


-0.03, 0.05 


Avinj ~ i sin(O.lt) 


R16 


1.0 


1.0 


0.001 


0.024 


50.0 


10.0 


3.0 


1.1 


-0.007, 0.007 


Auinj ~ i sin(0.3t) 


R17 


1.0 


3.0 


0.01 


0.22 


200.0 


30.0 


5.0 


1.1 


-0.0035, 0.013 


At)i„j ~ i(sin0.5t)^ 


R19 


1.0 


1.0 


0.05 


0.15 


50.0 


5.0 


2.0 


1.1 


-0.035, 0.05 


Auinj ~ i(sin0.2t)^ 



the acceleration process. The positive forces are strongly 
localized across the relection shocks. However, we don't 
see there effect in a positive jet rotation with in the jet 
beam. In comparison, the negative forces are distributed 
over a wide range between the internal shocks. As a 
result the (negative) rotation velocity of the jet beam in- 
creases along the propagation direction. The increase is 
step-wise with each of the inter-shock region adding to 
the overall rotation (see the step-wise increase of dark 
color along the jet flow in the upper right Fig. 111. 



Along the jet beam the negative Lorentz force wins and 
the overall jet rotation is negative. The situation is some- 
what different for the material which is entrained along 
the jet flow. Here the negative Lorentz force is more 
strongly localized and the positive force wins resulting 
in an overall positive rotation of the ambient (entrained) 
flow. 

This picture of a gradual increase of rotation by cu- 
mulative action of Lorentz forces along the path of the 
material is also supported by the small rotation feature 
visible near the jet head. A closer look shows that the re- 
spective Lorentz force and velocity patterns are slightly 
shifted. This tells us the jet ^-acceleration is not, say in 
situ, but that what we observe in rotation is the result of 
previously acting Lorentz forces. The same trend can be 
seen in the small scale toroidal velocity structures around 
the jet in the positive velocities where observe a similar 
shift between the positive Lorentz force and the positive 
toroidal velocity. 



7.4. Time evolution of angular momentum 



The magnetic fleld structure resulting from our rather 
basic and general simulation setup evolves into a complex 
helical fleld distribution implying a similarly complex 
electric current system. With the propagation of the out- 
flow and the evolving jet-internal dynamics the toroidal 
and poloidal field components change direction resulting 
in a similar variation of the torques along and across the 
jet. This is demonstrated exemplarily in Fig. l6j where 
we show the distribution of the specific toroidalLorentz 
forces, resp. the acceleration in toroidal direction. Com- 
parison of positive (right) and negative (left) Lorentz 
forces, highlights the filamentary structure of the force 
distribution. The structural details clearly depend on 
the specifics of the simulation setup and deserve a closer 
look by a future investigation. Observational indications 
for such a filamentary structure may be obtained from 
radial velocity profiles across the jet if the spatial reso- 
lution will be sufficient. 

The time-evolution shown in Fig. [6] is typical for the 
early jet evolution - the jet is penetrating the ambient 
medium and thereby changes the dynamical state com- 
pared to its initial injection profile. However, jet material 
and ambient gas are coupled by the magnetic field which 
is anchored into both components. The differential ro- 
tation between jet and ambient gas results in shear and 
torques. In later evolutionary stages when the jet has 
bored its funnel through the ambient gas, the differen- 
tial rotation decreases, and the twist in the field structure 
is reduced. Jet injections of different velocity will then 
lead to internal shocks (see simulation R16, Fig. [sjjo]). 

Figure 12 reveals the time sequence of the angular mo- 
mentum evolution as domain-integrated value for sim- 
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Figure 3. Axisymmetric 1.5D simulation of a propagating MHD 
shocis. Simulation RJ21 considers a rotating jet injected into an 



ambient medium with 



jet 



-0.2, but otherwise the same pa- 



rameters as RJ18 (see Fig.[2]l. Shown is the shock evolution at 
dynamical time t = 2 with a grid resolution of 1.000 cells / unit 
length. The bottom panel shows the specific total angular momen- 
tum pv^ (solid), with the magnetic (dashed), and hydrodynamic 
(dotted) contribution. 



ulation ROl. Shown arc the different contributions to 
the total angular momentum budget (i.e. not their spe- 
cific value). Since this is an density- weighted value, it 
is observationally more relevant than the specific value 
(simply speaking, more gas will radiate more). The ob- 
served jet rotation will not only depend on the rotation 
speed alone, but also on how much material is i n fact 
rotating with such speed. The example of Fig. 12 shows 
at time t = 20 the kinetic angular momentum on average 
more positive than negative. This suggests that one may 
observe on average a positive rotation of the jet although 
the negative velocity distribution is broader (see Fig. [gI) . 

At dynamical time t = 45 the jet bow shock leaves the 
computational domain. The jet material set into rotation 
by the processes we have discussed above will propagate 
out of the computational domain. Since simulation ROl 
has a time-independent injection boundary condition, no 
new shocks will be initiated which could enforce rotation 
and the overall jet rotation will fade away - as indicated 
by Fig. U 



As for simulation ROl we show in Fig. |13| the time evo- 
lution for the different contributions to tEe angular mo- 
mentum budget for simulation R12. In difference to ROl 
the amount of angular momentum is constantly growing. 
This is due to the presence of the steady reflection shock 
which continue to convert magnetic angular momentum 
to kinetic angular momentum, even when the bow shock 
has left the computational domain. 

8. CONCLUSIONS 

We have performed axisymmetric MHD simulations of 
jet propagation into an ambient gas applying the PLUTO 
code. The goal of our study was to investigate how the 
gas flow in jets and outflows is set in rotation by Lorentz 
forces arising MHD shocks. 

The process we discuss relies on a helical magnetic 
field structure which carries magnetic angular momen- 
tum from the jet base to the asymptotic jet. This could 
be converted into kinetic angular momentum by action of 
a toroidal torque which is induced by shock compression 
of the toroidal field. To demonstrate the applicability of 
the principal mechanism, for simplicity we have injected 
a non-rotating jet material into the ambient gas for most 
cases. 

Our simulations show that the proposed mechanism 
is indeed feasible, and accelerates the jet in toroidal di- 
rection by Lorentz forces which are induced by a shock 
compressed toroidal magnetic field. Depending on the 
magnetohydrodynamic parameters of the jet flow, the re- 
sulting rotational velocities are in the range of 0.1%— 1% 
of the jet propagation speed. 

Jets with a high external Mach number gain the high- 
est rotation speed. This is because the mechanism works 
most efficient in the termination shock. 

The resulting velocity field is highly complex result- 
ing in a filamentary structure in the toroidal torque and 
toroidal velocity distribution. The material which is set 
in rotation across the terminal shock, flows back respec- 
tive to the propagating jet and forms a cocoon with a 
highly tangled rotation pattern. 

Interestingly, of jets with the same initial magnetic 
angular momentum budget (i.e. same ~ iJ^B^), those 
with a low poloidal field (low magnetic flux) experience a 
stronger toroidal torque, and, thus gain higher rotational 
speed. This is a consequence of the higher compression 
rates due to the weaker poloidal field (more J-type than 
C-type shock). 

For the jets with steady injection, the toroidal velocity 
field decays as soon as the bow shock has left the com- 
putational domain. On the long term the injected gas 
flows along the channel bored by the jet during previous 
time steps and only weakly interacts with the ambient 
medium. For jets with time-dependent injection, further 
shocks are generated at a constant rate maintaining the 
jet rotation. 

Jets with low internal Mach number develop standing 
oblique reflection shocks which constantly (but stepwise) 
accelerate the jet material in toroidal direction. 

In summary, by demonstrating the feasibility of 
toroidal jet acceleration by MHD shocks in a helical mag- 
netic field, we propose the following implications. 

i) Jets of young stars are typically seen in radiation of 
shocked gas. Thus we expect this material to be affected 
by toroidal torques along its path from the jet launch- 
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Figure 4. Axisymmetric 1.5D simulations RJ22 (top) and RJ23 (bottom) considering jets with same magnetic angular momentum ^ B^Bz 
(and Poynting flux), but different strength of the field components, = AB^ (RJ22), B^ = 45,^ (RJ23), respectively (with otherwise the 
same parameters). Shown is the toroidal velocity (left), and kinetic angular momentum (right) at dynamical time t = A. 
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Figure 5. Axisymmetric 1.5D simulations RJ22 (top) and RJ23 
(bottom) considering jets with the same magnetic angular momen- 
tum ~ B^Bz (and Poynting flux), but different strength of the 
field components, B^ = AB^ (RJ22), B^ = AB^ (RJ23), respec- 
tively (with otherwise the same parameters) . Shown is the specific 
total angular momentum (solid), with the magnetic (dashed), and 
hydrodynamic (dotted) contribution. 

ing region close to the star up to asymptotic regime of 




a collimated jet. To extrapolate from observed angular 
velocities far from the source to the angular momentum 
of the jet material close to its origin is therefore ques- 
tionable. 

ii) Jet rotation, as observationally indicated in jets 
from young stars is feasible. Jets emerge from rotating 
disks and are ejected with kinetic angular momentum. 
In addition to this, MHD shocks could additionally ac- 
celerate the jet material by further converting magnetic 
to kinetic angular momentum. 

iii) Similar processes can be expected from relativis- 
tic stellar jets and extragalactic jets, leading to a rapid 
rotation in these jets. 

The results presented in this paper are derived from 
kinetic / dynamical modeling. Preliminary simulations 
considering cooling indicate not major effect, although 
we see that cooling results in a somewhat stronger shock 
compression and therefore slightly higher maximum ro- 
tation speeds. Detailed comparison with observations 
would require also modeling of radiation processes (radia- 
tive MHD), or post-processing of emission maps, which 
is beyond the scope of this paper. 
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APPENDIX 
COMPARISON SIMULATIONS 



Here we show comparison simulations applying a similar set of parameters (see T ab.|4t as 'Ryu & Jones ( 1995 1. Figure 
14 shows simulation RJll which basically repeats the simulation shown in Ryu fc Jones^ (1995[ ), l''ig. 2a, tor which we 
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Figure 6. Axisymmetric 2.5D simulation run ROl considers a Iiigh (external) Mach number jet with weak magnetic field. Sliown is (from 
top left to bottom right) the toroidal velocity at dynamical time t = 10,30, and for the latter also the axial velocity Vz, the toroidal field 
strength B^, and the negative and positive toroidal Lorentz force tFj^^^ (on a log scale). Shown is the inner part r < 5.0 of the whole grid 
of (r < 7.0,2; < 30.0). The inner, equidistant grid consists of (100 X 810) cells for (r < 1.5, z < 30.0). Attached to that is a logarithmically 
scaled grid for r > 1.5. The jet nozzle is located between < r < 1.0. 



have applied Cartesian coordinates as well. Figure 15 shows simulation RJ09 with similar dynamical parameters, but 
run in a cylindrical coordinate system, which can be considered more appropriate to our general aim of investigating 
jet rotation. In the latter case we have considered only one perpendicular vector component. The initial setup with 
vanishing radial magnetic field and vanishing radial motion is very well conserved. 
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Figure 11. Axisymmctric simulation run R12. Shown are the density distribution (top left), the radial velocity distribution (lower left), 
the negative (top) and positive (bottom) toroidal specific Lorentz force components (middle), and negative (top) and positive (bottom) 
toroidal velocity components (right) at dynamical time t = 50. 
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Figure 12. Time evolution of the domain-integrated angular mo- 
mentum (a.m.) contributions for example simulation ROl. Shown 
is the kinetic a.m. (positive contribution, thin dashed), kinetic a.m. 
(total contribution, thick dotted- dashed) , magnetic a.m. (negative 
contribution, dashed), kinetic a.m. (negative contribution, thin 
dotted- dashed), and the total a.m. {solid line). Around time t ~ 45 
the jet terminal shock leaves the computational domain. 
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Figure 13. Time evolution of the domain-integrated angular mo- 
mentum (a.m.) contributions for example simulation R12. Shown 
is the kinetic a.m. (positive contribution, thin dashed), kinetic a.m. 
(total contribution, thick dotted- dashed), magnetic a.m. (negative 
contribution, dashed), kinetic a.m. (negative contribution, thin 
dotted- dashed), and the total a.m. {solid line). 
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Figure 14. Simulation RJll in 1.5D Cartesian coordinates, equivalent to Fig. 2a of |Ryu &c Jones ]|T995]|. Shown are: density; gas pressure; 
axial and radial velocities; axial and radial magnetic field strength (from top left to bottom right] at dynamical time t = 10 (grid resolution 
5000 cells per physical length 2.0). 
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Figure 15. Simulation RJ09 in cylindrical coordinates and in axisymmetry, similar to Fig. 2a of |Ryu fc Jones] l |1995| |. Only one of the 
perpendicular vector components - the toroidal component - is non-vanishing. Shov^n are: density; gas pressure; axial, radial, and toroidal 
velocities; axial, radial and toroidal magnetic field strength (from top left to bottom right) at dynamical time t = 10 (grid resolution 3000 
cells per physical length 2.0). 



